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Abstract 

Stochastic network calculus provides an elegant way to characterize traffic and service processes. However, little effort has 
been made on applying it to multi-access communication systems such as 802.11. In this paper, we take the first step to apply 
, it to the backlog and delay analysis of an 802.11 wireless local network. In particular, we address the following questions: In 

' applying stochastic network calculus, under what situations can we derive stable backlog and delay bounds? How to derive the 

backlog and delay bounds of an 802.11 wireless node? And how tight are these bounds when compared with simulations? To 
answer these questions, we first derive the general stability condition of a wireless node (not restricted to 802.11). From this, 
■j^ , we give the specific stability condition of an 802.11 wireless node. Then we derive the backlog and delay bounds of an 802.11 

■ node based on an existing model of 802.11. We observe that the derived bounds are loose when compared with ns-2 simulations, 

' indicating that improvements are needed in the current version of stochastic network calculus. 
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I. Introduction 



Network calculus provides an elegant way to characterize traffic and service processes of network and communication 
systems. Unlike traditional queueing analysis in which one has to make strong assumptions on arrival or service processes 
(e.g., Poission arrival process, exponential service distribution, etc) so as to derive closed-form solutions yj, network calculus 
^ . allows general arrival and service processes. Instead of getting exact solutions, one can derive network delay and backlog 
bounds easily by network calculus. Deterministic network calculus was proposed in [2| [3| |4| I'Sl, etc. However, most traffic 
and service processes are stochastic and deterministic network calculus is often not applicable for them. Therefore, stochastic 
network calculus was proposed to deal with stochastic arrival and service processes Jl] Q Jll flU ifTol ifTTI lfT2l . 

There have been some applications of stochastic network calculus [TTl |14| [15| |T6'|. However, little effort has been made 
^*) ' _ on applying it to multi-access communication systems. In the paper, we take the first step to apply stochastic network calculus 
to an 802.11 wireless local network (WLAN). In particular, we address the following questions: 
\ • Under what situations can we derive stable backlog and delay bounds? 
Tij" ■ • How to derive the backlog and delay bounds of an 802.11 wireless node? 
\ • How tight are these bounds when compared with simulations? 

In this paper, we answer these questions and make the following contributions: 

• We derive the general stability condition of a wireless node based on the theorems of stochastic network calculus. From 
this, we give the specific stability condition of an 802.11 wireless node. 

■ • We derive the service curve of an 802.11 node based on an existing model of 802.11 IfTS I. From the service curve, we 
^ \ then derive the backlog and delay bounds of the node. 

• The derived bounds are loose in many cases when compared with ns-2 simulations. We discuss the reasons and point out 
future work. 

This paper is organized as follows. In Section [III we give a brief overview of stochastic network calculus. In Section |III1 
we present the stochastic network calculus model of a wireless node. In Section HVl we derive the general stability condition 
of a wireless node. In Section [V] we derive the backlog and delay bounds and the stability condition for an 802.11 node. 
In Section |Vl] we compare the derived bounds with simulation results. Related work is given in Section IVIII and finally. 
Section [Villi concludes the paper and points out future directions. 

II. Stochastic Network Calculus 

In this section, we first review basic terms of network calculus and then cite the results of stochastic network calculus which 
we will use in this paper There are various versions of arrival and service curves. We adopt virtual backlog centric (v.b.c) 
stochastic arrival curve and weak stochastic service curve in our analysis. 



A. Basic Terms of Network Calculus 



We consider a discrete time system where time is slotted (t — 0,1,2,...)- A process is a function of time t. By default, 
we use A{t) to denote the arrival process to a network element with A{0) = 0. A{t) is the total amount of traffic arrived to 
this network element up to time t. We use A*{t) to denote the departure process of the network element with yl*(0) ~ 0. 
A*{t) is the total amount of traffic departed from the network element up to time t. Let T (F) represents the set of non- 
negative wide-sense increasing (decreasing) functions. Clearly, A{t) G F and A* it) £ F. For any process, say A{t), we define 
A{s, t) = A{t) — A{s), for s < t. We define the backlog of the network element at time t by 

Bit)^A{t)-A*it), (1) 

and the delay of the network element at t by 

D{t) ^M{t : A{t) < A*{t + T)}. (2) 
Fig. [U illustrates an example of A{t) and A*(t) with B{t) and D{t) at t = 10. 
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Fig. 1. Illustration of A{t), A*{t), B(t) and D(t) 

In deterministic network calculus, A{t) can be upper-bounded by an arrival curve. That is, for all < s < t, we have 

A{s,t) < a{t - s), 

where a{t) is called the arrival curve of A{t). 

A busy period is a time period during which the backlog in the network element is always nonzero. For any busy period 
{to,t], suppose we have 

A*{t)-A*ito)>Pit-to), 

which means that the network element provides a guaranteed service lower-bounded by (3{t — to) during the busy period. We 
can let to be the beginning of the busy period, that is, the backlog at to is zero or A* {to) = A{to)- Therefore, 

A*{t)~A{to)>f3{t~to). 

The above equation infers A*{t) > info<s<f [A{s) + (3{t — s)], which can be written as 

A*{t)>A®l3{t), (3) 

where (g) is called the operator of min-plus convolution and fiit) is called the service curve of the network element. 

B. Results of Stochastic Network Calculus 

We cite the following definitions and theorems from jS] ifTTl except that we define Definition |6] by ourselves. 

Definition 1 (virtual-backlog-centric (v.b.c) Stochastic Arrival Curve): A flow is said to have a virtual-backlog-centric (v.b.c) 
stochastic arrival curve a <^ F with bounding function f ^ F, denoted by A ^vb< f,ct >, if for all t > and all x > 0, 
there holds 



P{ sup [A{s, t) - a{t - s)] > x} < f{x). 

0<s<t 



(4) 



Originally, in deterministic network calculus, we have A{s,t) < a{t — s) for all < s < t. However, there is usually some 
randomness in stochastic arrival processes and A{s,t) may not be upper-bounded by any arrival curve deterministically (e.g., 
traffic arrivals in {s,t] can be arbitrarily large in Poisson process). Thus, v.b.c stochastic arrival curve is proposed to tackle 
this problem. Roughly speaking, A{s,t) can exceed a{s,t) by x, but its probability is upper-bounded by f{x) which is a 
decreasing function of x. 

Deflnition 2 (Weak Stochastic Service Curve): A server S is said to provide a weak stochastic service curve /3 ^ J- with 
bounding function g E J", denoted by S ^ws< 9, /? >, if for all < > and all a: > 0, there holds 

P{A(g>P{t)~A*{t)>x}<g{x). (5) 

In deterministic network calculus, we have A*{t) > A(E) l3{t), which means that there is a service guarantee denoted by the 
service curve /3{t). However, there is usually some randomness in stochastic service process and thus a server may not always 
provide a guaranteed service curve deterministically. Thus, weak stochastic service curve is proposed to tackle this problem. 
Roughly speaking, A*{t) — A^ f3{t) can be less than ~x, but its probability is upper-bounded by g{x) which is a decreasing 
function of x. 

The utility of the above definitions is that if we can characterize the traffic by a v.b.c stochastic arrival curve and the server's 
service process by a weak stochastic service curve, then we can calculate backlog and delay bounds of the network element 
by Theorem [U 

Theorem 1 (Backlog and Delay Bounds): Consider a server fed with a flow A. If the server provides a weak stochastic 
service curve S ^ws< to the flow and the flow has a v.b.c stochastic arrival curve A ^vb< f, oi >, then 

(i) The backlog B{t) of the flow in the server at time t satisfies: for all i > and all a; > 0, 

P{B{t) >x}<f®g{x+ inf [P{s) - a{s)]). (6) 

s>0 

(ii) The delay D{t) of the flow in the server at time t satisfies: for ah i > and ah x > 0, 

P{D{t) >x}<f(g> 5(inf - ais - x)]). (7) 

By definition a{x) = when a: < in this theorem. Note that as noticed recently by researchers of network calculus, the 
formula of delay bound in this theorem often returns trivial results, which we will see in Section [V] 

By now, we have reviewed the key results of stochastic network calculus. Next, we will show how to calculate v.b.c stochastic 
arrival curve and weak stochastic service curve. 

In ifTTl . the author presented a theorem to facilitate calculation of stochastic arrival curves. Before showing the theorem, we 
first introduce {a{9), p{d))-upper constrained jS]. 

Definition 3 ({a{9), p{9))-upper constrained): A process A is said to be (ct(6'), p(0))-upper constrained (for some 9 > 0), 
if for all < s < we have 

^logEe'^^^^'^ <p{9){t-s) + a{9). (8) 

This definition is equivalent to Ee^"*'"^*^ < e^p(^)(*-^)+^'^(^), which means that A{s,tys moment generating function is 
upper-bounded. Two related concepts are defined as follows. 

Definition 4 (9-MER / 9-ER): A process ^'s minimum envelope rate (MER) with respect to 9 (9-MER), denoted by p*{9), 
is defined as follows: 

p* {9) = lim sup ^ sup log Ee''^(''"+*) . (9) 

t^co 9t s>o 

We say that A has an envelope rate (ER) with respect to 9 (0-ER), denoted by p{9), if p{9) > p*{9). 
The following theorem expresses the relationship between 9-ER. and ((t((?), /9(0))-upper constrained. 
Theorem 2 (Relationship of 9-ER and {(j{9), p{9))-upper constrained): 

(i) If A is (CT(6'),/9(6'))-upper constrained, then p{9) is 9-ER of A. 

(ii) If A has 9-ER. p{9) < oo, then for every e > there exists (t^{9) < cxo so that A is {(t^{9), p{9) + e)-upper constrained. 
Now we have two kinds of traffic characterization: v.b.c stochastic arrival curve and {<t{9), p{9))-\ippei constrained. The 

following theorem establishes the connection between them. 

Theorem 3 (v.b.c Stochastic Arrival Curve of {a{9), p{9))-upper constrained Process): Suppose A{t) is {a'{9), p{9))-\ippeT 



constrained, then it has a v.b.c stochastic arrival curv^H A ^vb< /, a >, where 



a{t) = r-t 

0^(9) 

for any r > p{9) and a; > 0. 

This theorem indicates that if we can show that the traffic is {a {9), p {6}) -upper constrained, then we can get its v.b.c 
stochastic arrival curve by Eq. ( fTOl l. 

We now introduce the concept of stochastic strict server. This concept was inspired by the observation that a wireless 
channel can be described by an ideal service process and an impairment process. As we will see in Section |III] a wireless 
node can be modeled as a stochastic strict server 

Definition 5 (Stochastic Strict Server): A server S is said to be a stochastic strict server providing stochastic strict service 
curve f3 E with impairment process / to a flow iff during any backlogged period (s, t], the output A*{s, t) of the flow from 
the server satisfies 

A*{s,t)>p{t-s)-I{s,t). (11) 

We can easily find the weak stochastic service curve of a stochastic strict server by the following theorem. 

Theorem 4 (Weak Stochastic Service Curve of Stochastic Strict Server): Consider a stochastic strict server S providing a 
stochastic strict service curve $ with an impairment process /. If the impairment process has a v.b.c stochastic arrival curve, 
or / ^vb< g, 7 >, and j3 E T, then the server provides a weak stochastic service curv^ S '^ws< 9, 1^ > with 

m^m^i(t). (12) 

So far, we have cited all results of stochastic network calculus which we will use in this paper Finally, we define stable 
backlog and stable delay. A natural definition is to check whether the expectation of backlog (or delay) is finite. 
Definition 6 (Stable Backlog/Delay): The backlog B{t) is stable, if 

EB{t) < oo, Vi. (13) 

Similarly, the delay D{t) is stable, if 

ED{t) < oo, Vt. (14) 
We say that the backlog (or delay) bound of stochastic network calculus is stable if they can derive stable backlog (or delay). 

III. Stochastic Networli Calculus Model of a Wireless Node 

In this section we model a wireless node (not restricted to 802.11) by stochastic network calculus. In general, we can define 
one slot (t = 1) to be any duration of time and measure traffic amount in any unit (e.g. bits, bytes or packets). 

We consider a wireless node. Let A{t) denote the traffic arrived at the node from the application layer Suppose A is 
(o'a(0i), PA (^i)) -upper constrained. From Theorem|3] we have A ^vb< f,ct >, where 



a{t) = VA-t 



for any ta > Pa{Oi). 

We can model a wireless node by a stochastic strict server. Let the channel capacity be c traffic unit per slot. The departure 
process A*{s,t) — ${s,t) — I{s,t) during any backlogged period {s,t], where = c ■ t is the ideal service curve and / 
is the impairment process due to backoff, channel sharing and transmission errors. Since I{s, t) < c ■ {t ~ s), I has a finite 
d-MER. From Theorem|2] there exist (7/(^2) and pi{92) so that / is (ct/ (6*2 ), p/(<?2)) -upper constrained. Based on Theorem|3] 
we have / ^^f,< 5, 7 >, where 

l{t) = ri-t 

'The original theorem (Theorem 5.1 II II ) established a similar connection between maximum-backlog-centric (m.b.c) stochastic arrival curve and 
((t(S), p(6))-upper constrained, which is wrong as noticed recently by researchers of network calculus. However, one can easily see the theorem holds 
for v.b.c stochastic arrival curve. 

-In the original theorem (Lemma 4.2 111]), the impairment process has a m.b.c stochastic arrival curve. However, one can easily see that the theorem also 
holds for the impairment process which has a v.b.c stochastic anival curve. 



for any r/ > pi{92). 

From TheoremlH the node provides a weak stochastic service curve S ^ius< g,P >, where 

m = {c-ri)-t, (17) 

for any c > rj. 

Furthermore, from Theorem [T] we must have a{t) < fi{t), or equivalently, 

rA<c^ri. (18) 
Thus, P{B{t) > x} < f®g(x). Otherwise, if a{t) > I3{t), we get a trivial backlog bound, P{B{t) > a;} < f®g{-oQ) = oo. 

IV. Stability Condition of a Wireless Node 

One fundamental question we need to address is under what condition we can get stable B{t) and D{t) from stochastic 
network calculus, i.e., EB{t) < oo and ED{t) < oo. Before presenting our result, we first define the concept of envelop 
average rate. 

Definition 7 (Envelop Average Rate): The envelop average rate of a process A, denoted by a^, is defined as 

EA{s,s + t) 

ttA = lim sup . (19) 

Let a A and a/ be the envelop average rate of A and /, respectively. The following proposition shows the stability condition. 
Proposition 1 (Stability Condition): A wireless node has stable backlog and stable delay if 

aA < c~ aj. (20) 

Proof: We have shown that P{B{t) > x} < f ® g{x) if Eq. ( fTsT i holds. Thus, for any t. 



EB{t) = ^P{B(^)==^ + l}•(^ + l) 

i=0 

oo 

< ^P{i?(t)>z}.(z + 1) 

Z^O 
OO 

< 5I/'^5«-(* + l)<oo. (21) 

i=0 

Since f{x) and g{x) are exponentially decreasing functions according to Eq. ( fTSl l and Eq. ( fT&t , f ®g{x) is an exponentially 
decreasing function. Thus we have, for any t, EB{t) < oo. It is easy to see that for any t, ED{t) < oo. Otherwise, the service 
time is oo and thus limt_j.oo Ei?(t) = oo which contradicts Eq. jTH . 

Now we examine Eq. (fTSl l. From Eq. (flST l and Eq. (fT6] l. va — pa{Oi) + e and r/ — pj(^92) + e for any e > 0. Thus, Eq. (fTSl l 
is equivalent to 

PA{0i)<c^pi{92)~2e. 

From Theorem|2] we have pa{9i) — P*a{9i) + ei and pi{02) — P*j{S2) + ei for any ei > 0, where p\{Oi) and p*j{02) are 
6'-MERs of A and /, respectively. Equivalently, we have 

p*Ai0i)<c-p}{92)-2{e + ei). (22) 

Using Taylor expansion on ^^(6*1), we have 

P*a{0i) = lim sup -5-suplogEe^i'^(''''+*) 

= lim sup - — sup logE{l + 9iA{s,s + t) + 0{9l)) 

t^oo Oit s>o 

= lim sup -5- sup log (1 + eiEA{s, s + t) + 0{9l)) 
= lim sup -5- sup [6'iEA(s, s + i) + 0(6*?)]. 



Therefore, 



lim p\{9i) — lim sup 

01— !-0 t^-OO ^.>Q 



EA{s, s + 1) 



t 



= OA- 



Similarly, 



lim p*j{92) = aj. 



Therefore, there exist 6i and 62 so that p\(Oi) < + £2 and p*j{d2) < a/ + £2 for any enough small 62 > 0. So Eq. ( l22b 
is satisfied if 



Remarks: Since the proof is based on theorems of stochastic network calculus, it indicates that we can get stable backlog and 
delay bounds by stochastic network calculus as long as the condition of Eq. ( l20l ) holds. Since this condition is very general, 
we conclude that theoretically stochastic network calculus is effective. 



In this section, we apply the results in the previous section to calculate backlog and delay bounds for an 802.11 WLAN 
node. For simplicity, we assume there are n identical stations (or nodes) sending packets to an access point. All nodes operate 
in Distributed Coordination Function (DCF) mode with RTS/CTS turned off ifTTl . We consider an ideal channel, that is, 
transmission errors are only caused by collisions. Two packets are collided if their transmissions overlap in time. Besides, we 
assume that all DATA packets are of the same size. 

A. 802.11 DCF 

A node with a DATA packet (or simply packet) to transmit first monitors the channel activity. If the channel is idle for a 
period of time equal to a distributed interframe space (DIES), the node transmits. Otherwise, if the channel is sensed busy 
(either immediately or during the DIES), the node backs off, in which the node defers channel access by a random number of 
idle slots within a contention window (CW), ranging from to CW — 1. When the backoff counter reaches zero and expires, 
the node can access the channel. During the backoff period, if the node detects a busy channel, it freezes the backoff counter 
and the backoff process is resumed once the channel is idle for a duration of DIES. To avoid channel capture, a node must wait 
a random backoff time between two consecutive new packet transmissions, even if the channel is sensed idle for a duration of 
DIES. Once the packet is received successfully, the receiver will return an ACK after a short interframe space (SIES). Note 
that SIES is shorter than an idle slot so that there is no collision caused by a DATA packet and an ACK. 

802. 1 1 uses the truncated exponential backoff technique to set its CW. For example, in 802.1 lb, the initial CW is CWmin = 
32. Each time a collision occurs, CW doubles its size, up to a maximum of CWmax = 1024. When the packet is successfully 
transmitted, CW is reset to CWmin- The packet is dropped when it is retransmitted for six times and still not transmitted 
successfully. Eig.|2] shows some parameters of 802.11b used in our paper. 



QA < c - ai - 2{£ + ei + €2). 



Since e, ei and £2 can be arbitrarily small, the above equation is satisfied when 



aA < c — ai. 



V. 802.11 Backlog and Delay Bounds 
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Fig. 2. 802.11b parameters 



The duration of an ACK is the duration of a PHY header and an ACK header transmitted at basic rate, i.e., '"'^^^^a^'^ = 304/is. 
The duration of a DATA packet is the duration of an PHY header transmitted at basic rate plus the duration of an MAC header 
and its upper-layer payload transmitted at data rate. For example, suppose the upper-layer payload is 256 bytes, then the 
duration of the DATA packet is ^ + ^^nao«^'^ = 398.5Ais. 

B. Service Curve 

Since we only consider packets of equal size, we can measure traffic amount in packets. For simplicity, we measure time 
duration (e.g. SIFS, DIFS, DATA and ACK) in idle slots and we define that one slot in network calculus {t = 1) is equal to 
L idle slots, where 

L = DIFS + DATA + SIFS + ACK. (23) 

Note that sometimes "idle slot" only refers to a time period and it may not be idle. To avoid confusion in the following context, 
we will use "idle slot" (italic) to denote that the "idle slot" is indeed idle. 

In practice, it is difficult to calculate the impairment process / accurately since / depends on the complex interactions of 
traffic arrival and DCF. In this section, we perform the worst case analysis based on an existing model of 802.11 [18]. 

We assume that the system is in saturated condition, that is, the backlog at each node is always nonzero. Let r denote the 
transmission attempt probability per idle slot by a node and let 7 denote the conditional collision probability given that there 
is a transmission. We assume 7 is constant and independent for each transmission. Intuitively, this assumption becomes more 
accurate when the number of nodes n increases. In ifTSl , the authors derived a general formula relating r to 7, which is 

r= l + 7 + 7^ + -+7^ (24) 
^0 + 7&1 + 7^^2 + ••• + l^b^ 

This equation can be explained as follows. The numerator is the expected number of transmission attempts of a packet. The 
denominator is the expected total backoff duration (in idle slots) of a packet, where bi is the mean backoff duration plus 1 
(the 1 refers to the first idle slot of packet transmission) after the zth collision. In 802.11, bi — ^ •C'^""" where < i < 6. A 
packet suffering six collisions will be dropped from its buffer. In our model, we do not consider packet drops, which yields 
upper bounds of backlog and delay. 

The independence assumption of 7 implies that each transmission "sees" the system at steady state. Therefore, each node 
transmits with the same probability t. This yields 

7 = l-(l-Tr-i. (25) 

Combining Eq. (l24l l and Eq. dZSl ), we can solve t and 7. 

We introduce the following notations for an 802.11 node. The probability of no transmissions at an idle slot, denoted by 
Pnt, is (1 — r)". The probabiUty of having at least one transmission at an idle slot, denoted by Pt, is 1 — Pnt- The probability 
of a given node starting a successful transmission at an idle slot, denoted by Ps, is r(l — 7). For the given node, the probability 
of the other nodes starting transmissions at an idle slot, denoted by Po, is Pt — Ps — 7. 

Fig. [3] plots Eq. (l24l) in dashed line and Eq. dZSl l in solid line when n = 10, 20 and 100. The points of intersection are the 
solution of 7 and t for different n. When n increases, 7 increases and t decreases. Consequently, Pnt and Po increases, but 
Ps decreases. Therefore, the assumption of saturated condition gives the worst-case analysis. 

An 802.11 node can be seen as a stochastic strict server Clearly, the stochastic strict service curve /3(t) = t, which means 
that one packet is transmitted during one slot in the ideal case. In order to characterize the impairment process /, it is crucial 
to know Ee^^'^'''''"'''^ which we calculate as follows. 

We consider the transmissions of a given 802.1 1 node during t slots. From Eq. ( |23] ). there are tL idle slots in t slots, indexed 
1, 2, tL. At the first or the last idle slot, the transmission (if any) can be incomplete, that is, the transmission can start before 
the first idle slot or it can end after the last idle slot. We assume that at the first slot there is always a complete transmission 
by any of the other nodes except the given node, which actually overestimates Ee^^'-*'^"'"*''. Suppose there are i complete 
transmissions and zero or one incomplete transmission within the remaining {t — 1)L idle slots. The incomplete transmission 
occupies the last k idle slots where 0<fc<L— 1(A: = means that the last transmission is actually complete). Thus, there 
are {t — i — 1)L — k idle slots of no transmissions within the t slots. The probability of i complete transmissions and the 
last k idle slots occupied by an incomplete transmission, denoted by pk,i, is C(t-i-i)L-k+i-^nt ' Pt - Furthermore, the 

probability of the given node having j successful complete transmissions on condition that there are i complete transmissions, 
denoted by j, is C/ (Ps/ PtY (Po/PtY^-' ■ Finally, we have Ee^^^''''*+*'' is upper-bounded by 

■ E E Ep'^-^P..-'^*"^'^ + E Epo,.P,,,e«('-^). (26) 

fe=l i=0 j=0 i=0 j=0 
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Fig. 3. The Plots of Eq. (24) and Eq. (25) 

The first term is for the case that the last transmission is incomplete and the second term is for the case that the last transmission 
is complete. 

In general, we do not have an analytical form of Eq. (|26] |. so we resort to numerical methods and use Algorithm 1 to obtain 
ai{d) and pi{0) (see Appendix A-1). This algorithm is immediately inspired from Definition [3] Then we can use Eq. ( fT6l l and 
(ITtI i to obtain the node's weak stochastic service curve. 

Scenario 1: 

10 identical nodes sending packets to one access point 
The payload of a DATA packet is 256 bytes 

Fig. 4. Parameters of Scenario 1 



We illustrate the above calculations in Scenario 1 (see Fig. [Hi. From Eq. (l24l i and dZSl l, r = 0.037 and 7 = 0.293. Thus, 
Pnt = 0.680, Pt = 0.320, P, = 0.027 and Po = 0.293. Let M{t) be the value of Eq. Fig. |5] shows M(t) and its 

(j{0), {p{9))-nppei constrained curve when 9 — 1. We have p{l) — 0.948 and cr(l) = 0.096, as calculated by Algorithm 1. 

Furthermore, from Eq. (fT&t and ( fTTj i we have S ^ws< 9-, P > where 

p0.096 

9{x) 



I _ g0.948-r/ ' 

for any 1 > r/ > 0.948. Fig.|6]plots and g{t) when r/ = 0.968. 

C. StabUity Condition 

Before calculating the bounds, we first derive the stability condition. Since f3{t) = t, the channel capacity c = 1 packet 
per slot. Because 1 — a/ is actually the percentage of the node's successful transmission time, by Proposition [H the stability 
condition of 802.11 is 

D. Backlog and Delay Bounds 

We can immediately calculate backlog bounds by applying Eq. (flSTl-lfTTli into Theorem [T] The only technical issue is to 
select proper 9i and 92 a obtain tight bounds. Clearly, f{x) depends on 6i and r^, and g{x) depends on 62 and r/. According 
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Fig. 5. Scenario 1: the (o"(l), p(l))-upper constrained curve of a node's M{t) 
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Fig. 6. Scenario 1: a node's Weak stochastic service curve 

to Eq. (fTSb . we have ta < l — rj and P{B{t) > x} < f ®g{x). In addition, (ri) should be set as large as possible because 
f{x) {g{x)) decreases with ta [rj). Considering the above conditions, we have 

P{B{t) > x} < niin[/ g) g{x)] 
subject to 

rA > PAiOi),ri > pi{92) and + rj = 1. (28) 

In general, we do not have an analytical solution of min[/ ® g{x)] and we resort to numerical methods and use Algorithm 2 
to get a near-optimal solution (see Appendix A-2). 



As noticed recently by researchers of network calculus, the delay bound in Theorem [T] often returns trivial results. In our 
model in Section Hill it is easy to see that P{D{t) > x} < / g(0) = 1. We propose the following way to estimate delay 
bound. Little's law states that the average number of customers in a queueing system is equal to the average arrival rate of 
customers to that system, times the average time spent in that system ||T|. Let the average arrival rate is A. Assume the system 
can reach steady state when t ^ oo. Then we have the average backlog is limt_j.oo EB{t) and the average delay of each packet 
is greater than or equal to limj^oo ED{t) (by its definition in Eq. (|2]i, D{t) can be less than the delay of the bottom-of-line 
packet at t). Therefore, by Little's law, we have 

limt_^oo EB{t) 



lim ED{t) < 



A 



Finally, we apply Markov's inequality to the above equation and we have 

\imt^ooEB{t) 



P{ lim D{t) >x}< 



(29) 



(30) 



Besides, according to Eq. dlB, EB{t) < X^^o P{B{t) >«}■(« + 1). And we can use Eq. (|28]l to bound P{B{t) > i}. Note 
that Eq. ( l30t is derived when t -> oo. In practice, we can use this result to estimate delay bound when t is sufficiently large. 

VI. Performance Evaluation 

In this section, we use ns-2 simulations to verify our derived backlog and delay bounds for Poisson and constant bit rate 
(CBR) traffic arrivals. We carry out all experiments for Scenario 1 (Fig. Each simulation duration is 100 seconds (s) 
which is long enough to let a node transmit thousands of packets. Each data point (e.g. P{B{t) > x}) is calculated over 100 
independent simulations. 

A. Poisson Traffic 

Let A be the average traffic rate (packets/slot). In this case, we have ga — ^ (see Definition |7]i. For Poisson traffic, we have 



oo 

s.s+t) ^ ^ 

! = 



(At)' 



(31) 



where -^^^e is the probability of i packets arriving within (s, s + t]. From the above equation, Poisson traffic is (0, 



He"-!) 



upper constrained and we can obtain its arrival curve by Eq. dTsT l. 

From Eq. dZTl i. backlogs are stable when A < 0.079 (packet per slot). In Fig. |7] we plot the average backlog E[i3(i)] at 
t — 50s and A — 0.077, 0.079 and 0.081 in ns-2 simulations. We observe that there is sudden jump when A = 0.081, indicating 
the critical point of stability is indeed around 0.079. This figure also indicates the accuracy of the 802.11 model in Eq. (|24|) 
and Eq. 
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Fig. 7. EB(t) (t = 50s) when A = 0.077, 0.079, 0.081 



Experiment l(Scenario 1 with low Poisson traffic load) We set A = 0.04 to simulate low traffic load. 

Fig. [Hla) shows P{B{t) > x} in ns-2 simulations and Fig. |8lb) shows the upper bound of P{B{t) > x} calculated by 
Eq. ( |28] |. Note that stochastic network calculus gives very loose upper bounds. There may be two reasons. One is that we use 
the worst-case analysis in deriving the weak stochastic service curve of 802.11. The other is that there are many relaxations 
in proving the theorems of stochastic network calculus fTT]. For example, relaxations are used in deriving f{x) in Theorem[3] 
IT2J, and this theorem is popularly used in deriving arrival curves and service curves (see Section HIH). The first reason may 
not be the key reason because we will see in Experiment 2 the bound is even looser when we increase arrival rate and 
make the channel near saturated. The second reason seems to be the key reason. We will see in Experiment 3 that backlog 
bounds improve substantially for CBR traffic where we are able to derive the arrival curve by hand without using Theorem [3] 
This indicates that refinements are needed in stochastic network calculus so as to tighten the bounds. Moreover, we found 
that backlog bounds are sensitive to adjusting parameters (i.e., di, 62, ta and r/). So it is necessary to use Algorithm 2 to 
minimize the bounds. 

We also conduct simulations to verify delay bounds at i = 50.s. Since the backlog bounds are too loose, in order to avoid 
trivial validation, we use E_B(t) in ns-2 simulations to validate Eq. (|29] l and Eq. dSOl l (assume that t = 50s is sufficiently large 
so that we can apply these equations). Actually, we get ED{t) < 0.0205s by Eq. ( |29] |, which tightly bounds ED{t) = 0.0186s 
in ns-2 simulations. Fig. [Sic) shows P{D{t) > x} in ns-2 simulations and Fig. [Ud) shows the upper bound of P{D{t) > x} 
calculated by Eq. (l30l l. Clearly, P{D{t) > x} is upper-bounded by Eq. 
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Fig. 8. Experiment 1: When t = 50s (a) P{B(t) > x} (b) upper bound of P{B{t) > x} (c) P{D{t) > x} (d) upper bound of P{D(t) > x} 



Experiment 2(Scenario 1 with high Poisson traffic load) We set A — 0.07 to simulate high traffic load. 

Fig. |9fa) shows P{B{t) > x} in ns-2 simulations and Fig. |9jb) shows the upper bound of P{B{t) > x} calculated by 
Eq. ( |28] l. In this case, 9i (62) is much smaller than that of Experiment 1 so as to satisfy the constraint in Eq. ( |28] l, resulting 
in looser f{x) and g{x). Therefore, stochastic network calculus gives further loose backlog bounds. 

We also conduct simulations to verify delay bounds at t = 50s. Since the backlog bounds are too loose, in order to avoid 
trivial validation, we use EB{t) in ns-2 simulations to validate Eq. (|29] l and Eq. ( |30] | (assume that t — 50s is sufficiently large 
so that we can apply these equations). Actually, we get ED{t) < 0.0299s by Eq. (|29] l, which tightly bounds ED{t) = 0.0296s 
in ns-2 simulations. Fig. |9lc) shows P{D{t) > x} in ns-2 simulations and Fig. ^d) shows the upper bound of P{D{t) > x} 
calculated by Eq. (|30l i. Clearly, P{D{t) > x} is upper-bounded by Eq. (|30] i. 



(a) (b) 




0.02 0.04 0.06 0.08 0.1 ^ 0.02 0.04 0.06 0.08 0.1 

X (second) x (second) 



Fig. 9. Experiment 2: When t = 50s (a) P{B{t) > x} (b) upper bound of P{B{t) > x] (c) P{D(t) > x} (d) upper bound of P{D{t) > x} 



B. CBR Traffic 

Let A be the average traffic rate (packets/slot). In this case, we have ga — A. It is easy to see that supo<5<t [A{s, t) — X ■ 
{t — s)] < 1 for all t because packets arrive one by one in a constant time interval. Thus, we have A ^vb< f,ct> where 

a{t) = Xt 

m - {?; ll\ (32) 

Again, from Eq. (IZTT i. the stability condition is A < 0.079 (packet per slot). In Fig. [TOl we plot the average backlog ¥\B{t)] 
al t = 50.S for A — 0.077, 0.079 and 0.081 in ns-2 simulations. We observe that there is sudden jump when A = 0.081, 
indicating the critical point of stability is indeed around 0.079. This figure also indicates the accuracy of the 802.11 model in 
Eq. (EllianddBll. 

Experiment 3(Scenario 1 with low CBR traffic load) We set A — 0.04 to simulate low traffic load. 

Fig. [TlT a) shows P{B{t) > x} in ns-2 simulations and Fig. [TTT b) shows upper bound of P{B{t) > x} calculated by 
Eq. (l28T l. The backlog bounds are much tighter in CBR traffic than those in Poisson traffic (see Experiment 1). The main 
reason is that we can derive a tight f{x) by hand instead of by Theorem [3] 

We also conduct simulations to verify delay bounds at t = 50s. Since the backlog bounds are still loose, in order to avoid 
trivial validation, we use EB{t) in ns-2 simulations to validate Eq. (|29] | and Eq. ( l30l l (assume that t = 50s is sufficiently large 
so that we can apply these equations). Actually, we get ED{t) < 0.0090s by Eq. (|29] |. which tightly bounds ED{t) ~ 0.0089s 
in ns-2 simulations. Fig. fTTT c) shows P{D{t) > x} in ns-2 simulations and Fig. fTTT d) shows the upper bound of P{D{t) > x} 
calculated by Eq. (l30l l. Clearly, P{D{t) > x} is upper-bounded by Eq. (|30] |. 

Experiment 4(Scenario 1 with high CBR traffic load) We set A — 0.07 to simulate high traffic load. 

Fig. fTSl a) shows P{B{t) > x} in ns-2 simulations and Fig. fT2l" b) shows the upper bound of P{B{t) > x} calculated by 
Eq. ( |28] |. The backlog bounds are much tighter in CBR traffic than those in Poisson traffic (see Experiment 2) because f{x) 
is tight here. 

We also conduct simulations to verify delay bounds at t = 50s. Since the backlog bounds are still loose, in order to avoid 
trivial validation, we use EB{t) in ns-2 simulations to validate Eq. ( |29] l and Eq. dSOl l (assume that t = 50s is sufficiently large 
so that we can apply these equations). Actually, we get ED{t) < 0.0274s by Eq. ( |29] |. which is close to (although does not 
bound) ED{t) = 0.0296s in ns-2 simulations. Fig. \l2\ c) shows P{D{t) > a;} in ns-2 simulations and Fig. [T2l d) shows the 
upper bound of P{D{t) > x} calculated by Eq. ( l30l l. Clearly, P{D{t) > x} is upper-bounded by Eq. (l30l l. 

To sum up, the current version of stochastic network calculus often derives loose bounds when compared with simulations, 
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Fig. 10. EB{t) {t = 50s) when A = 0.077, 0.079, 0.081 
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Fig. 11. Experiment 3: Wlien t = 50s (a) P{B{t) > x} (b) upper bound of P{B{t) > x} (c) P{D(t) > x} (d) upper bound of P{D{t) > x} 

especially in the case of high traffic load. Therefore, stochastic network calculus may not be effective in practice. 

VII. Related Work 

In this section, we first present relate work on stochastic network calculus and then on the performance analysis of 802.11. 

The increasing demand on transmitting multimedia and other real time applications over the Internet has motivated the study 
of quality of service guarantees. Towards it, stochastic network calculus, the probabilistic version of the deterministic network 
calculus |T| |3| [Tl (5], has been recognized by researchers as a promising step. During its development, traffic-amount-centric 
(t.a.c) stochastic arrival curve is proposed in |6|, virtual-backlog-centric (v.b.c) stochastic arrival curve is proposed in [9J and 
maximum-backlog-centric (m.b.c) stochastic arrival curve is proposed in ITTl . Weak stochastic service curve is proposed in Q 
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Fig. 12. Experiment 4: When t = 50s (a) P{B(t) > x} (b) upper bound of P{B{t) > x} (c) P{D{t) > x} (d) upper bound of P{D{t) > x} 



||8l and stochastic service curve is proposed in flO\. In flT], Jiang showed that only the combination of m.b.c stochastic arrival 
curve and stochastic service curve has all five basic properties required by a network calculus (i.e., superposition, concatenation, 
output characterization, per-flow service, service guarantees) and the other combinations only have parts of these properties. 
Jiang also proposed the concept of stochastic strict server to facilitate calculation of stochastic service curve. Moreover, he 
presented independent case analysis to obtain tighter performance bounds for the case that flows and servers are independent. 
However, there are a few bugs in his results recently found by researchers of network calculus, such as the trivial delay bound 
in Theorem 3.5 111] and Theorem 5.1. Therefore, we adopt v.b.c stochastic arrival curve and weak stochastic service curve in 
our study since we only consider backlog and delay bounds (i.e., service guarantee), ignoring the other properties. 

There have been some applications of stochastic network calculus. In |13|, Jiang et al. analyzed a dynamic priority 
measurement-based admission control (MBAC) scheme based on stochastic network calculus. In |14|, Liu et al. applied 
stochastic network calculus to studying the conformance deterioration problem in networks with service level agreements. In 
1151, based on stochastic network calculus, X. Yu et al. developed several upper bounds on the queue length distribution of 
Generalized Processor Sharing (GPS) scheduling discipline with long range dependent (LRD) traffic. They also extended the 
GPS results to a packet-based GPS (PGPS) system. Finally, Agharebparast et al. modeled the behavior of a single wireless 
link using stochastic network calculus [16|. However, little effort has been made on applying stochastic network calculus to 
multi-access communication systems such as 802.11. 

Existing work on the performance of 802.11 has focused primarily on its throughput and capacity. In |fT9l , Bianchi proposed 
a Markov chain throughput model of 802.11. In [18], Kumar et al. proposed a probability throughput model which is simpler 
than Bianchi's model. In our paper, we adopt Kumar's model to derive the service curve of 802.11. There are also some work 
on queueing analysis of 802.11. In ll20l . Zhai et al. assumed Poisson traffic arrival and proposed an M/G/1 queueing model 
of 802.11. More generally, Tickoo proposed a G/G/1 queueing model of 802.11 f 2T\ ||22| . To our best knowledge, we are the 
first to model the queueing process of 802.11 based on stochastic network calculus. 

VIII. Conclusion and Future Work 

In this paper, we presented a stochastic network calculus model of 802.1 1. From stochastic network calculus, we first derived 
the general stability condition of a wireless node. Then we derived the stochastic service curve and the specific stability condition 
of an 802.1 1 node based on an existing model of 802.1 1. Thus, we obtained the backlog and delay bounds of the node by using 
the corresponding theorem of stochastic network calculus. Finally, we carried out ns-2 simulations to verify these bounds. 

There are some open problems for future work. First, we derived the service curve based on an existing 802.11 model. Thus, 
the accuracy of the service curve depends on the accuracy of the model. An open question may be whether we can derive 
the service curve of 802.11 without using any existing models. Second, we assumed the worst-case condition (i.e., saturate 



condition) in our analysis. Can we remove this conservative assumption? Besides, under the worst-case assumption, we can 
assume flows and servers are independent and perform independent case analysis obtaining tighter backlog and delay bounds. 
This is also one of our future work. Third, we observe that the derived bounds are loose when compared with ns-2 simulations, 
calling for further improvements in the current version of stochastic network calculus. 
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Appendix A-1: Algorithm 1 (Numerical Calculation of cfi{0) and pi{0)) 

1) Let M{t) ~ sup3>g{i logEe^^(*^*+*)}. Obviously, M{t) is an increasing function of t with A/(0) ~ 0. We define axes 
t and axes t±_ (vertical to t) on a plane, and plot M (t) on it. We define the slope of M (t), s{t) — M (t) — AI {t — 1). 

2) We calculate s(t) for i = 1, 2, 3, ... until it converges at some t*, i.e., (1 - e) ■ s{t* - 1) < s{t*) < (1 + e) • s(t* - 1) 
where e is a small number, e.g. 10~^. 

3) We draw a straight line l{t) with the slope s{t*) crossing the point (t* , M{t*)) . Obviously, the line crosses the point 
(O, M{t*)~s{t*)t*y The maximum displacement between M{t) and l{t) (in the direction of t±), Vm — niaxo<t<t* {M{t) — 
l{t)}. We shift l{t) by v„i in the direction of t± and get l{t). Clearly, l{t) upperbounds M{t). In other words, we have 
pj{9) = s{t*) and ai{9) = M{t*) - s{t*)t* + v^- 

Appendix A-2: Algorithm 2 (Numerical Calculation of Eq. ^28^ ) 

In each iteration, we generate a sample of Oi, 62, ta and r/. If they satisfy the condition of Eq. ( |28] ), we calculate 
min[/ (g) g{x)] for the current and past iterations until it converges. Sample generations can use the interpolation or Monte 
Carlo method over valid ranges of the variables. 



